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We show that autocatalytic networks of e-machines and their population dynamics differ substan- 
tially between spatial (geographically distributed) and nonspatial (panmixia) populations. Gener- 
ally, regions of spacetime-invariant autocatalytic networks — or domains — emerge in geographically 
distributed populations. These are separated by functional membranes of complementary e-machines 
that actively translate between the domains and are responsible for their growth and stability. We 
analyze both spatial and nonspatial populations, determining the algebraic properties of the auto- 
catalytic networks that allow for space to affect the dynamics and so generate autocatalytic domains 
and membranes. In addition, we analyze populations of intermediate spatial architecture, delineat- 
ing the thresholds at which spatial memory (information storage) begins to determine the character 
of the emergent auto-catalytic organization. 
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I. INTRODUCTION 

Almost certainly, the sophisticated mechanisms for 
self-replication found today in living cells were not 
present in the earliest replicators [THS]. Instead, some 
hypothesized that the first replicators were autocatalytic 
networks of functional molecules that collectively were 
capable of self- reproduction [4j [5]. Numerous studies 
have been devoted to these and analogous networks with 
the hope of understanding pre-biotic evolution [SHE]- A 
particular class of such models represented network el- 
ements with mathematical constructs, including regular 
expressions, the A-calculus [5], and, recently, e-machines 
[TU] , e-Machines are especially useful since they support 
well defined and computable measures of structural com- 
plexity pi] and, equally important, these measures ex- 
tend directly to networks of interacting e-machines. 

Early models often assumed a Turing gas in which ev- 
ery network clement has an equal chance of interacting 
with every other, thereby ignoring a network's spatial 
configuration. However, natural and engineered evolu- 
tionary systems are not architected this way since ele- 
ments of physical systems always have some spatial re- 
lationship that determines which components interact. 
This observation led to studies of spatial pattern forma- 
tion in evolutionary and autocatalytic systems [T2HH]- 

Here, we consider networks of single-state 
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e-machines — a simplifying initial focus that helps 
to highlight the role of population architecture. We 
show that the behavior of spatially distributed popula- 
tions differs substantially from that of the nonspatial 
Turing gas. More to the point, we determine the 
algebraic properties of the network that lead to the 
emergence of distinctive organizations. 

II. BACKGROUND 

e-Machines serve as models of stochastic finite and 
infinite computation |15) . An e-machine consists of 
a set S of causal states and a set of transitions be- 
tween those states. In the particular variation used here 
(e-transducers), each transition is labeled by both an in- 
put symbol x £ A and an output symbol y G A. An 
e-machine in state a £ S reads an input symbol x and 
chooses from all transitions from a the one labeled x. 
The e-machine follows the chosen transition to a next 
state a' € S while emitting the output symbol y corre- 
sponding to that transition. In this way, e-machines can 
be viewed as mapping an input language to an output 
language, perhaps probabilistically. 

Following Ref . [10] , we focus on single-state e-machincs 
over the input-output alphabet A = {0, 1}. (Results for 
multi-state e-machines appear in a sequel. Our goal here 
is to highlight the effects of space, uncomplicated by the 
richness that comes with using multi-state e-machines.) 
There are 16 such e-machines and we denote the set of 
all them by T. Each e-machine can be represented as a 
2x2 binary matrix M, where = 1 means that the 
e-machine reads in symbol i — 1 while emitting symbol j — 
1. We number the 16 e-machines — T = {T , . . . ,Ti 5 } — 
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by finding the decimal equivalent of the binary number 
AI11M12M21M22 for each e-machine. Thus, for example, 
e-machine Tn has the matrix representation: 



and set: 



M = 



1 
1 1 



(1) 



e-Machine pairs interact in a population by compo- 
sition, under which T is closed and forms a monoid 
[16j . With the matrix representation, e-machine 
composition — Tj, o T a — is simply matrix multiplication 
where, after multiplying, any positive matrix element is 
set to 1. From this, it is straightforward to compute the 
interaction matrix M, where T a o Tt, = T c if and only if 
M Q+ i i6+ i = c, where a, b, c € {0, ... , 15}: 
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III. e-MACHINE SOUPS 

A population, or simply a soup, T is a configura- 
tion of an n x n regular toroidal lattice. At each time 
t = 0, 1,2, . . ., every lattice location contains a sin- 
gle e-machine, denoted r* ■ 6 T. The population size is 
N = n 2 . Each location (i,j) is initialized to contain an 
e-machine uniformly chosen at random from T. 

The population dynamics is specified by e-machinc 
composition, interaction, and update. We define a func- 
tion 0:TxTxT^Tby: 



6(T a ,T b ,T c ) = 




if T a oT c ^ T 
otherwise 



(3) 



We write 9{A,B,C) for sets A,B,C C T as shorthand 
for the set {6(a,b, c) : a € A,b € B,c € C}. We must 
use 9 to prevent the non-e-machine To from being pro- 
duced since the fact that T o T a = T a o T = T for ah 
T a e T implies that if T could be produced, it comes to 
dominate any population. 

Each time step t we choose a location at random 




with probability 
with probability 
with probability 
with probability 



(4) 



and set T^ 1 = T\ d for all (k, I) ^ (i,j). Thus, two verti- 
cal or horizontal neighbors to Tij are chosen, composed, 
and the e-machine resulting from their composition is 
used to replace r* if it is not T . This replacement 
scheme is meant to be locally analogous to the replace- 
ment scheme used in Ref. [10], which will facilitate direct 
comparisons in the following. 

In addition, at each time step a certain amount of dif- 
fusion of e-machines occurs due to spatial mixing. For 
this, let be a Gaussian distribution with variance v 
and mean 0. At each time step, c e-machines are cho- 
sen at random in T and each is swapped with a random 
e-machine at a distance chosen from the distribution £„. 
For more generality, we allow c to be any real number 
and swap c e-machines per time step on average. 

When v and c are large, there is considerable spatial 
mixing and one expects the dynamics to behave like a 
nonspatial population, where e-machines have an equal 
chance of interacting with every other. However, when v 
and c are small, there is little spatial mixing and, as we 
shall see, the population dynamics change substantially. 

To summarize, as a stochastic dynamical system the 
soup's state at time t is the population's configuration 
r*. While Eq. Q determines the local probabilistic up- 
date at each site, we use 6 to formally denote the global 
(probabilistic) update for the entire configuration over 
one time step: 



r* +1 = e o r* 



(5) 



And so, one goal is to understand the trajectories 
{r°, r 1 , . . .}. Another is to analyze the structure inside 
the r*. For this, we use or to denote a spatial shift of 
the soup configuration: 



T' = a ft oT 



(6) 



where n = (Ai, Aj) is the vector by which the configu- 
ration is shifted horizontally and vertically: 



= r 



(i+Ai) mod n,(j+A_7") mod n 



(7) 



Finally, we view the population either as a configu- 
ration r of spatially positioned e-machines or as a dis- 
tribution f of e-machine types without regard to spatial 
location. The fractions f = (f 1 , . . . , f 15 ) of e-machines of 
type T a on T are P(r) = Pr(T Q e T). 
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IV. THE PANMIXIA SOUP 

The least-spatial architecture of T, where c « N and 
v fa n, was previously studied in Ref. [10]. This is the 
case of a panmixia population, in which all e-machines 
can interact with any other. When the size of T is large, 
f(i) = f(r*) follows a simple population dynamics. 

First, for large N the discrete-time dynamic is well ap- 
proximated by continuous time, since each discrete time 
updates only a single lattice location and this means 
that changes to T* and so to i(t) are relatively small 
(oc A^ 1 ). Second, the probability of adding an e-machine 
Ti is that of picking two e-machines T a and T b , such that 
T a ° Tf, — Ti, times the probability that the e-machine 
replaced is not T,. Third, the probability of removing an 
e-machine T t is that of picking two e-machines, T a and 
T b , such that T a o T b ^ T t and also T a o T b ^ T , times 
the probability of picking a Ti to replace. Thus, the rate 
of change of f l (t) is given by the differential equation: 



dP 
~dt 



= i-r 



{ a { b - f i 

T a oT b =Ti 



J2 f a f b 



(8) 



T a oT b ^Ti 
T a oT b =jtT 



where the sums run over all pairs satisfying their sub- 
scripted condition. That is, ^ T ot 6 =t runs over all or- 
dered pairs (T a ,T b ) such that T a °o T b = T % . Equation ([8} 
also determines the steady-state probability distributions 
of e-machine types, which are simply solutions of: 



rff 

dt 



= . 



(9) 
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FIG. 1: Panmixia population evolution: e-Machine-type dis- 
tribution f as a function of replication time (t/N), with 
n = 10 3 , N = 10 6 , c = 3000, and v = 1000. 



Figure [T] shows the distribution f(t) over e-machine 
types as a function of time (replication) for a simulation 
with c = 3000 and v — 1000. This is essentially the 
same population behavior reported in Ref. |10j . It was 
confirmed numerically that Eq. ^ closely predicts the 
results in Fig. [l] In Fig. [TJ all 16 single-state e-machines 
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T3,T 5 ,T W , T12 
Ti,T 2 ,T4, T g 

rri rri rri ryi rri rr\ 

j 6j j 9, J 11, J 13, J 14 


Fast Growth 
Medium Growth 
No Growth 
Fast Decay 



TABLE I: e-Machine-type behavior classes in the panmixia 
soup of Fig. [T] 



are shown. They partition themselves into four classes 
such that those in the same class behave similarly; see 
Table [I] Six e-machines die away, while nine persist in 
a closed, self-maintaining, and dynamically stable meta- 
machine, as shown in Ref. [10] . 



V. GENERAL REPLICATORS 

While these e-machine classes can be understood by 
examining the solutions of Eq. ^ as done in Ref. [10] , 
the results can also be directly predicted by examining 
the algebraic structure of the set T under the composi- 
tion operator determined by M. And this observation is 
critical to predicting the behavior of alternate population 
architectures. First, we address the panmixia population. 
The central idea is to find subsets of e-machine types that 
map onto themselves under the population dynamics; in 
other words, to identify M-invariant subsets of T. 

Definition 1. A set S of e-machines is a general repli- 
cator ( GR ) if for all a G T we have: 



and 



9(T,a,S)U0(S,a,T) CSU{a} 



0(T,a,S)U0{S,a,T) ^ {a} 



This parallels the definition of an ideal in semigroup 
theory, except that we must be more careful (and less 
elegant) with the definition since 9 is not a binary relation 

In T there are many GRs, including T itself. We in- 
tentionally used 9 to exclude the non-e-machine T from 
being considered a GR. By design, To cannot be pro- 
duced. GRs are important because they are the simplest 
form of replicator. Suppose S is a GR. As evolution 
progresses, elements in the soup are replaced with some 
element of 9(T, a, S) or 9(S, a, T). By the definition of a 
GR, though, 9(T, a, S) and 9(S, a, T) contain elements of 
S and so generally we expect elements of S to produce 
more elements of S. 

However, not all GRs are equal since some are subsets 
of others and therefore are produced more readily. There 
are, in fact, many GRs in T, but there is one GR, call 
it i7, that is minimal in that no subset of Q is a GR, 
but every GR contains il. (Note that, generally, such a 
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set is not guaranteed to exist in a semigroup with a zero 
element [16].) 

Proposition 1. fi = {T u T 2 , T 3 , T 4 , T 5 , T 8 , T 10 , T 12 , T 15 } 
is t/ie minimal general replicator in T. 

Proof. By direct verification of Def. [7] and observing 
that fi, short any one of its e-machines, is not a GR. 

In the panmixia population, we expect that elements 
of 51 come to dominate the soup. And this is what is 
observed in the simulations (Fig. [I] and Table |l| . While 
Eq. (JsJ) explains the specific values of f , an understanding 
of M's structure leads to a direct explanation for why the 
set Fast Decay is removed from the soup: None of the 
e-machines in Fast Decay are in Q. 

VI. THE SPATIAL SOUP 

The population that we consider next is that of a spa- 
tially configured soup where c = or v = 0. In these 
limits, no spatial mixing is added to the system and the 
model behaves essentially as an asynchronous, probabilis- 
tic cellular automaton. Specifically, each element is re- 
placed by the composition of two of its neighbors if that 
composition does not result in T . The pair of neighbors 
composed is chosen with uniform probability, according 
to Eq. @. 



0.45 | r 

0.4 - 




5 10 15 20 25 30 35 40 45 50 



tlN 

FIG. 2: Spatial population dynamics: e-Machine-type dis- 
tribution f as a function of replication time with n = 10 3 , 
N = 10 6 , and c = 0. 

Figure [2] shows the e-machine-type distribution f (t) for 
c = and N = 10 6 . Note that initially the population 
behaves quite similarly to the panmixia case. This is to 
be expected since the simulations begin with identical 
initial configurations so at t = the range of possible 
interactions is effectively the same for both architectures. 
As in the panmixia population, elements not in Q are 
quickly removed from the soup. 

By t/N ps 3.4, however, the populations start to be- 
have differently For example, the e-machinc T15 that 
is most readily made in the panmixia case is also read- 
ily made in the spatial case, until t/N — 3.4, when it 
begins to be removed from the population. Again, the 



e-Machine Type 


Early 


Late 


Ti r 

1 15 

T3,T 5 , 7j.O) T\2 
T 2 ,T 4 
Ti,T 8 

f~J~l rri rri rri rj~l rri 

1 &, ±7, + 9, J 11 ) J 13, J 14 


JJdolj \jlUW Lill 

Medium Growth 
Logistic Growth 
No Growth 
Fast Decay 


vj\jK,c% y 

Decay 
Saturation 
No Growth 
Removed 



TABLE II: e-Machine-type behavior classes in the two- 
dimensional spatial population for early and late times for 
Fig. [2] Early times: t/N < 3.4; late times: t/N > 3.4. 



e-machines partition themselves into subsets whose ele- 
ments behave similarly. At early times before t/N = 3.4, 
these are identical to those seen in the panmixia pop- 
ulation except that T 2 and T4 undergo logistic growth, 
while Ti and Tg do not. Table |TT| summarizes the overall 
behavior, for both early and late times. 

More striking than the evolution of the e-machine-type 
distribution f(£), however, are the spatial patterns that 
emerge in T* . Figure [3] shows a 70 x 70 region of T* at 
increasing times and one 500 x 500 region at a late time. 
Regions of T 2 and T4 form stable domains that grow to 
dominate the soup. Within a region of T 2 or T4 there 
exist e-machines that cannot be replaced, typically due 
to an "elastic" collision — an interaction producing the 
non-e-machine To that, by 6, simply leaves the e-machine 
at the lattice location alone. 

When a region of T 2 meets a region of T4, e-machines 
Ti and Ts form on the boundary, since T 2 o T4 = Ti and 
T40T2 = T s . Moreover, TxoTx = 7\ and T 8 oT 8 = T s and 
so the boundaries are self-sustaining along the interface. 

The spatial soup exhibits many of the nontrivial spon- 
taneous patternings common to reaction-diffusion sys- 
tems that exhibit Turing instability [I7H2D] . Finding 
a set of reaction-diffusion partial differential equations 
equivalent to this model and the spatial analog of Eq. [8j 
however, remains an open problem. 

Over a large number of time steps, the spatial popu- 
lation looks essentially like that at t/N = 20.0 in Fig. 
[3] If enough time passes, though, the soup eventually 
divides itself into two connected regions of T 2 and of T 4 , 
or one will take over completely. This requires an ex- 
tremely large number of replications. Which e-machine 
region dominates in the long run appears to be randomly 
determined. The overall process is highly reminiscent of 
spinodal decomposition in which a mixed solution sepa- 
rates into stable component phases [21]. A more direct 
connection to the predictions of that theory awaits fur- 
ther effort. 

VII. SPATIAL-REPLICATOR DOMAINS AND 
THEIR MEMBRANES 

After time t/N pa 3.4, the characters of the spatial 
and panmixia populations begin to diverge substantially, 
with patterns emerging in the spatial T. Those patterns 
consist of domains of nearly homogeneous e-machines of 
type T 2 or of type T4. They are separated by domain 
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t/N = 9.0 t/N = 20.0 : 70 x 70 t/N = 20.0 : 500 x 500 



FIG. 3: Emergent spatial replicators and their membranes: A 70 x 70 region of T with n = 10 3 , N = 10 6 , and c = v = 0. The 
last image, however, is of a 500 x 500 region. Each color corresponds to one of the 16 e-machine types: T2 is green, T4 is dark 
gray, Ti is yellow, Tg is purple, and T15 is blue. 



walls that consist of e-machines T\ and T 8 . As we shall 
see, these walls play the role of functional membranes 
that actively translate between the domain e-machines. 

To summarize, then, the first observation about the 
spatial population is that it is completely differently orga- 
nized. The persistent set of e-machines differs markedly 
from those found with panmixia, as do their roles and 
interactions. Why did they emerge when the population 
is embedded in space? We can explain the emergence of 
these structures by defining a mechanism in the inter- 
action network M that affects the population dynamics 
when there is sufficient "spatial memory". This mecha- 
nism is called a spatial replicator (SR). 

Definition 2. A set S of e-machines is a spatial repli- 
cator (SR) if 

9{S,S,d(S, S, T)) C S, 
6(S, S,0(T,S, S)) C S, 
8{6(S, S,T),S, S) C S, and 
0(0(T,S,S),S, S) C S. 



Notice that this is similar to the definition of a GR, 
except that SRs require two applications of 9 to produce 
an element of S. Due to this, SRs only replicate when 
there is spatial information storage: If S is a SR, an 
element of S must compose with an e-machine in T to 
produce an e-machine T y . T y must then compose with 
an element of S again to produce a new element of S. 

In the panmixia population there is no notion of adja- 
cency. Therefore, T y does not have an increased chance of 
interacting with an element of S again. However, space 
allows a way for SRs to replicate: If by chance a suf- 
ficiently large domain of T consists only of elements of 
some SR S, e-machines on the domain border will be re- 
placed by elements of 6(S, S, T) or 6(T, 5, S). Later, these 
boundary e-machines will be replaced with elements of 

6(S, S,6(S, S, T)) U 6(S, S, 0(T, S, S)) 

U 6(6(S, S, T), S, S) U 0(0(T, S, S),S, S) , (10) 

which are all elements of S. Thus, a domain of S will 
grow. This argument lays out how the algebraic struc- 
ture that space adds to M leads to domain growth. In 
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FIG. 4: The effective e-machine interactions at various times. The nodes are the 15 possible e-machines. The darkness of their 
interiors indicates their relative fractions f % in the populations. The edges connecting them indicate the frequency of their 
interactions; darker is more frequent. 



addition, since 6(S, S, S) — S, domains of S are self- 
maintaining. 

As with a GR, one expects that if S' is a SR and is 
also a subset of some other SR S, then S' eventually 
will come to dominate the soup. For the spatial popu- 
lation of single-state e-machines, it is easy to check that 
{Ti, T 2 , T4, T 8 } is a SR. Dynamically, the membranes of 
T\ and T$ grow their respective domains of T4 and T-i- 
Unsurprisingly, then, T2 and T4 come to dominate T. In 
contrast, X15 — an e-machine made so readily by the pan- 
mixia population — begins to dominate the simulation un- 
til "spatial memory" begins to act. Beyond that point 
T15 is replaced by SRs. 

Figure [4] shows how these interactions develop over 
time by giving the effective e-machinc interactions on an 
undirected graph, where the darkness of the e-machinc 
nodes indicates their relative frequencies f 1 and the dark- 
ness of the connecting edges indicates the relative fre- 
quency of interactions. As a complement to this, Fig. 
[5] gives the e-machine interaction network at long times. 
This should be compared to the nine-e-machine meta- 
machine in the panmixia soup: see Fig. 3 of Ref. [10] . 

The explanation for why a spatial population organizes 




FIG. 5: The e-machine interaction network at long times in 
the spatial soup. An arrow labeled Tt going from node T a to 
node T c denotes T c = T b oT a . 



differently and, in effect, selects different individuals for 
different roles relies on SRs, a local mechanism. The do- 
mains and walls, however, are nonlocal structures. Given 
a local dynamic that produces them, how do we describe 
their global structure? 
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To do this we adapt the computational mechanics 
analysis of emergent domains and particles in deter- 
ministic, synchronous cellular automata (CAs) [52] to 
the stochastic, asynchronous spatial population dynam- 
ics here. Speaking informally, a domain is a patch of 
spacetime that has the same "texture" when shifted in 
time, space, or both. How much one must shift each do- 
main so that its texture maps onto itself is one crude 
measure of the structural complexity of the domain. 

Definition 3. A domain is a set of configurations T s — 
{Tjj G S}, where S is an SR, that satisfies: 

1. Temporal- shift invariance: T s = $T S , and 

2. Spatial-shift invariance: V s = a n T s for some spa- 
tial offset n. 

It is straightforward to see that, in the spatial popula- 
tion, homogeneous regions of SRs — {T 2 } and {T4} — are 
domains. Determining which additional e-machines — the 
nonreplaceable ones mentioned above — can be embedded 
in these domains, such that the regions are still domains, 
is a more difficult calculation that we will not attempt 
here. 

The membranes separating the domains are comple- 
mentary e-machine types that actively function to trans- 
late between domains on either side, being composable 
with the e-machine types in those domains. And, over 
time, sets of membrane e-machines map back into them- 
selves. In the spatial population, as noted above (Fig. 
|, Ti = T 2 o T 4 and T 8 = T 4 o T 2 . The definition of a 
membrane replicator and a membrane express these for- 
mally. 

Definition 4. A set M of e-machines is a membrane 
replicator (MR) if there are two domain SRs, S and S' , 
such that 

6(S,M,6(S,M,S')) C M, 

6(S,M,9(S\M, S)) C M, 
O(0(S, M, S'), M, S') C M, and 
6(6(S',M,S),M,S') C M. 



Notice that this is similar to the definition of an SR, 
except that MRs require the bounding domain SRs to 
produce e-machines in the MR. 

Definition 5. A membrane between two domain SRs S 
and S' is a set of configurations: 

T ss' = r s MT s' ^ ( U ) 

where M is a membrane replicator, that is temporal- shift 
invariant: 

r ss' _ $y ss ' . (12) 



Observe that membranes, like domains, are temporally 
invariant, but, unlike domains, they are not spatially 
shift invariant. This is what it means, in fact, for a struc- 
ture to have a location in space. 

Membranes translate between the, possibly distinct, 
kinds of information stored in the bounding SR domains. 
This view parallels that of, for example, the particles 
which perform computations in evolved cellular automata 
|23j . Specifically, when the spatial population dynamic 
is deterministic, then SR domains and membranes are 
directly analogous to the domains and particles of cel- 
lular automata computational mechanics [22 j. The par- 
ticular class of population dynamical systems considered 
here differ, however, in important ways: They have an 
asynchronous, stochastic spatial update. A more detailed 
framing of these differences and the methods to analyze 
them will be reported elsewhere. 

In summary, Defs. [3j |4j and [5] specify algebraic con- 
straints that, in principle, can be solved to find SR do- 
mains and membranes. Even in the simpler case of the 
computational mechanics analysis of deterministic CAs, 
solving the analogous set of constraints, though well de- 
fined, is at present difficult. And so, we will leave detail- 
ing the algorithms for these calculations to a sequel. 

VIII. INTERMEDIATE SPATIAL 
POPULATIONS 

The e-machine soups here were specifically designed 
to permit investigation of populations that are neither 
completely spatial nor completely nonspatial by vary- 
ing the mixing parameters c and v. To avoid present- 
ing extremely high-dimensional data, we pick a metric 
that distinguishes the spatial versus nonspatial popula- 
tion dynamics and observe how that metric varies with 
c and v. For the metric we use f 15 (£), the percentage of 
e-machines of type T15 in T*. In the nonspatial simula- 
tion, f 15 w 0.223 after t/N = 50 steps, while f 15 » 0.0 
after t/N — 50 in the spatial model. In this way, Fig. [6] 
shows the transition between spatial and nonspatial be- 
havior for the model as a function of c and v. 




FIG. 6: f (t) after t/N = 50 time steps as a function of c 
and v with N = 10 6 . 



The two plots in Fig. [7] provide a clearer demarcation 



8 



0.2 - ..... 

0.15 - 

0.1 - 

0.05 - 

Q 1 . . . . L . A J . J .l . l . lj 1 . . .. L . . 1 . -I- 1 l l l 1 1 I I I I I 

0.001 0.01 0.1 1 10 

V 

0.25 I i i q i i q i i q i i q i i q i i q r~ 

0.2 - 
0.15 - 

m 
cm 

0.1 - 

0.05 - 

L... J j...lJ. ...*•••'* I I 1 L_lJ 1 L_lJ 1 L_lJ 1 — 

0.001 0.01 0.1 1 10 100 1000 

c 

FIG. 7: Top panel: f 15 (i) after t/N = 50 time steps with 
N = 10 6 and c = 10.0. Bottom panel: f 15 after t/N = 50 
time steps with N = 10 6 and v = 1.0. Note that in both, the 
horizontal axis is logarithmic. 

of the transition at which spatial memory is lost by show- 
ing slices along the v and c axes. Notably, the plots show 
that the "effect of space" is logistic in logc or log v. This 
indicates that populations are sensitive to both parame- 
ters in similar ways, even though the physical realizations 
of varying c (number of e-machines swapped) and v (spa- 
tial range of mixing) are quite different. 

IX. CONCLUSION 



in the spatial case. Structural analysis of the interac- 
tion network, focusing specifically on its algebraic prop- 
erties, was essential to understanding the mechanisms 
that drive the spontaneous emergence of organization. 
To this end, we drew a useful, if preliminary, connection 
to the computational mechanics analysis of domains and 
walls structures in spatially extended systems. 

The membranes observed (and denned) here are not 
merely concentration gradients, as seen in familiar pat- 
tern formation systems. Rather, they are entities that, 
due to their interaction specificity, actively translate be- 
tween replicators in neighboring domains. The resulting 
spatial organization suggests that one gets spontaneous 
compartmcntalization without designing it in at the be- 
ginning. Recall that compartmentalization is often cited 
as one of the key early evolutionary steps on the road to 
increasing biological complexity [2]. It is not such a dif- 
ficult step, after all. The e-machine soup is a very simple 
system, built with a minimal set of physical, chemical, 
and biological assumptions. 

It is perhaps no surprise if we mention that analogous 
investigations of one-dimension e-machine soups leads to 
similar results. This and the two-dimensional, inter- 
mediate dimensional, and effectively infinite dimensional 
(panmixia) soups lead one to wonder about the role of 
the three-dimensions in which we know that life arose. Is 
there something special? 

Constructively, in concert with natural spontaneous 
pattern formation, evolution may very well commandeer 
structures of different spatial dimension as ways that in- 
sure various kinds of functionality or intrinsic computa- 
tion. We appreciate that evolution is opportunistic and 
takes advantage of nature's propensity to spontaneously 
organize. And so, is three-dimensional space more advan- 
tageous in, say, the range of structures that it supports, 
compared to one and two dimensions? Or does evolution 
take advantage of all possible dimensions? The empir- 
ical evidence speaks rather clearly, if not to the former 
question, then to the latter: genetic information is stored 
in one-dimensional structures, biological cells sport two- 
dimensional membranes, and organisms move and behave 
in three dimensions. 



The analysis demonstrated that the capacity for using 
spatial memory can be determined by examining the in- 
teraction network. Moreover, it also showed that the be- 
havior of populations which are neither completely spa- 
tial nor completely nonspatial exhibit a systematic de- 
pendence on the degree of mixing. One can imagine 
several ways to generalize the definition of SRs — as the 
definition presented above is sufficient, but by no means 
necessary — for the spatial populations to behave differ- 
ently from the nonspatial populations. 

While the exact population dynamics of f may be cap- 
tured by an ODE such as Eq. ([8| , this coarse-graining is 
not useful to understanding the population dynamics' key 
structural properties — the pattern formation highlighted 
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